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Field-dependent specific heat and neutron scattering measurements were used to explore the 
antiferromagnetic S=^ chain compound CuCb • 2((CD3)2SO). At zero field the system acquires 
magnetic long-range order below Tat = 0.93 K with an ordered moment of 0.44/i_g. An external field 
along the b-axis strengthens the zero-field magnetic order, while fields along the a- and c-axes lead 
to a collapse of the exchange stabilized order at fioH c = 6 T and fioHc = 3.5 T, respectively (for T = 
0.65 K) and the formation of an energy gap in the excitation spectrum. We relate the field-induced 
gap to the presence of a staggered g-tensor and Dzyaloshinskii-Moriya interactions, which lead to 
effective staggered fields for magnetic fields applied along the a- and c-axes. Competition between 
anisotropy, inter-chain interactions and staggered fields leads to a succession of three phases as a 
function of field applied along the c-axis. For fields greater than /j.qH c , we find a magnetic structure 
that reflects the symmetry of the staggered fields. The critical exponent, /3, of the temperature driven 
phase transitions are indistinguishable from those of the three-dimensional Heisenberg magnet, while 
measurements for transitions driven by quantum fluctuations produce larger values of f3. 

PACS numbers: 75.25.+Z, 75.10.Pq, 74.72.-h 



I. INTRODUCTION 

Weakly-coupled antiferromagnetic (AF) S — h chains 
can be close to one or several quantum critical points and 
serve as interesting model systems in which to explore 
strongly-correlated quantum order. An isolated AF S = 
\ Heisenberg chain, described by 

H = J 2_j ' Sn+i , (1) 

n 

is quantum critical at zero temperature. Due to the ab- 
sence of a length scale at a critical point, the spin corre- 
lations decay as a power law, and the fundamental exci- 
tations are fractional spin excitations called spinons that 
carry S = Interactions that break the symmetry of 
the ground state can induce transitions to ground states 
of distinctly different symmetry, such as one-dimensional 
long-range order at zero temperature in the presence 
of an Ising anisotropy or conventional three-dimensional 
long-range order at finite temperature in the presence of 
weak non-frustrating interchain interactions. 3 

Magnetic fields break spin rotational symmetry and 
allow tuning the spin Hamiltonian in a controlled way. 
When it couples identically to all sites in a spin chain, an 
external field uniformly magnetizes the chain and leads to 
novel gapless excitations that are incommensurate with 
the lattice. Experimentally this is observed with meth- 
ods that couple directly to the spin correlation functions, 
such as neutron scattering. Inelastic neutron scattering 
experiments on model AF S = h chain systems such as 
copper pyrazine dinitrate^ show that upon application of 
a uniform magnetic field there are gapless excitations at 



an incommensurate wave vector that moves across the 
Brillouin zone with increasing field.— This is expected 
to be a rather general feature of a partially magnetized 
quantum spin system without long range spin order at 
T = K. When even and odd sites of the spin chain 
experience different effective fields, there is an entirely 
different ground state which is separated by a finite gap 
from excited states. This has been observed in magnetic 
materials with staggered crystal field environments and 
Dzyaloshinskii-Moriya (DM) interactions. Most notably, 
the generation of an excitation energy gap was first ob- 
served in copper benzoate using neutron scattering^ and 
the entire excitation spectrum of the S = ^ chain in the 
presence of an effective staggered field was mapped out 
in CuCl 2 • 2DMSO (CDC)2. These experiments showed 
that the excitation spectrum in high staggered fields con- 
sists of solitons and breathers that can be described as 
bound spinon states, and that the one-dimensional char- 
acter of the spinon potential leads to high-energy exci- 
tations that are not present in the absence of staggered 
fields^ 

In this paper, we explore the strongly anisotropic field 
dependence of long range magnetic order in CDC. Us- 
ing neutron scattering, we have directly measured the 
order parameters associated with the various competing 
phases. We show that a magnetic field destroys the low- 
field long-range order that is induced by super-exchange 
interactions involving Cl~ cations, and that an applied 
field induces a first-order spin-flop transition followed 
by a second-order quantum phase transition. Field- 
dependent specific heat measurements provide evidence 
for a gapless excitation spectrum in the low-field phase 
and an energy gap for larger fields. These results com- 



bined with neutron measurements enable a comprehen- 
sive characterization of the dominant spin interactions ir 
this system, including a quantitative determination of tht 
staggered gyromagnetic tensor and the DM interactions 
which give rise to the field-induced gap. 

II. CRYSTAL STRUCTURE 

CDC crystallizes in the orthorhombic space grouf 
Pnma (No. 62) with room-temperature lattice constant; 
a = 8.054A, b = 11.546 A, and c = 11.367 A £ Basec 
on the temperature dependence of the magnetic suscep- 
tibility, CDC was proposed as an AF S=h chain witl 
exchange constant J = 1.43 meVJ^ Inspection of th( 
crystal structure (the positions of the Cu 2+ sublattic( 
are given in Appendix[AJ suggests strong super-exchangt 
interactions between Cu 2+ ions via Cl~ cations, so thai 
the Cu 2+ ions form spin chains along the a-axis, as showr 
in Fig. [TJ Inelastic neutron scattering measurements de- 
scribed below yield a consistent value of J = 1.46 meV. 

The local symmetry axis of the Cl~ environment lies 
entirely in the ac-plane, but alternates its orientation 
along the a-axis. This alternation of the local crystalline 
environment leads to an alternating g tensor in CDC of 
the following form: 




2.31 (-l)V \ 
2.02 =g u + (-ir g s , 

(-l)V 2.14 J 



(2) 

where n denotes the n-th spin in the chains along the 
a-directioni^. g u is the uniform and g s is the staggered 
part of the g factor. 

The crystal symmetry also allows the presence of DM 
interactions^^ 

H D M=^D„.(S„_ 1 xS n ). (3) 

n 

The DM vector D„ = (— l) n Dh alternates from one 
chain site to the next as can be seen from the space 
group symmetry: First, the structure is invariant under 
a translation along the a-axis by two Cu sites, meaning 
that there can be at most two distinct DM vectors. Sec- 
ond, the ac plane is a mirror plane, implying that the 
DM vectors point along the b-axis. Third, the crystal 
structure is invariant under the combined operation of 
translation by one Cu site along the chain, and reflection 
in the ab-plane. 

It was shown that alternating DM interactions can re- 
sult in an effective staggered field H st upon application 
of a uniform external field H fRef. ll3T). Taking into ac- 
count also the staggered g tensor, the effective staggered 
field that is generated by a uniform field can be written 
as 

g u H st = i-Dxg u H + g s H. (4) 




FIG. 1: Position of the Cu 2+ (red), CI" (green) and CT 
(blue) in the crystal structure of CDC. The dominant superex- 
change interactions are via Cu 2+ -Cl~- Cu 2+ paths forming 
linear chains along the a-axis. The axes of the quasi-planar 
CUCI2O2 groups alternates by 22° along the chain direction, 
leading to a staggered gyromagnetic tensor g. 



III. EXPERIMENTAL TECHNIQUES 

Single crystals of CDC were obtained by slow cool- 
ing from T = 50°C to 20°C of saturated methanol so- 
lutions of anhydrous CuCl2 and dimethyl sulphoxide in 
a 1:2 molar ratio. Emerald-green crystals grow as large 
tabular plates with well developed (001) faces. Specific 
heat measurements were performed on small, protonated 
crystals with a typical mass of 15 mg in a dilution refrig- 
erator, using relaxation calorimetry in magnetic fields up 
to 9 T applied along the three principal crystallographic 
directions. Measurements were performed in six differ- 
ent magnetic fields applied along the a- and b-axes, and 
eight different fields applied along the c-axis. 

Neutron diffraction measurements were performed on 
deuterated single crystals of mass lg using the SPINS 
spectrometer at the NIST Center for Neutron Research. 
Inelastic neutron scattering experiments were performed 
on multiple co-aligned single crystals with a total mass of 
7.76 g using the DCS spectrometer at NIST. The SPINS 
spectrometer was configured with horizontal beam colli- 
mations 45'/fti ( 58 Ni guide) -80' -80' -300'. A pyrolytic 
graphite (PG) monochromator was used to select inci- 
dent neutron energies of either Ei = 3 meV, 5 meV (us- 
ing the (0, 0, 2) reflection) or 20 meV (using the (0, 0, 4) 
reflection). For Ei = 3 meV and 5 meV a cooled Be 
filter was used before the sample to eliminate contam- 
ination from higher-order monochromator Bragg reflec- 
tions. The DCS measurements were performed with an 
incident energy £^=4.64 meV and an angle of 60° be- 
tween the incident beam and the a-axis, as detailed in 
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FIG. 2: Specific heat of CDC as a function of temperature for 
magnetic field applied along the a, b and c-axis. The peaks 
in the temperature dependence correspond to the onset of 
long-range magnetic order. The solid (and the mostly hidden 
dashed) line corresponds to the low-temperature specific heat 
calculated in the framework of the sine-Gordon (boson) model 
and described in Eq. [10] (Eq. [6}. Dotted line plotted in (a) 
corresponds to aT 3 lattice contribution to the specific heat. 



Ref. 0. 



IV. SPECIFIC HEAT MEASUREMENTS 

The temperature dependence of the specific heat at 
low fields has a well-defined peak for all field direc- 
tions, as shown in Fig. [2j The peak position varies with 
field strength and direction, all of which suggests that 
the peak corresponds to the onset of long-range mag- 
netic order. At zero-field, the transition temperature is 
T/v=0.93 K, in agreement with prior zero-field heat ca- 
pacity measurements ^ When the field is applied in the 
mirror plane, either along the a or the c-axis, the Neel 
temperature Tjv decreases with increasing field, as shown 
in Fig. [UJa) and Fig. [2c). The peak becomes unobserv- 
able above noH c =5.7 and /ioH c =3.9 T for fields applied 
along the a- and c-axis, respectively. At higher fields, 
the specific heat is exponentially activated, with spectral 
weight shifting to higher temperatures with increasing 
field. This is evidence for a spin gap for magnetic fields 
along the a- and c-axes. In contrast, a magnetic field 
along the b-axis enhances the magnetic order and the 




FIG. 3: Semilogarithmic plot of the specific heat C(H,T) 
in the high-field gapped phase as a function of the inverse 
temperature l/T for different magnetic fields applied along 
the (a) a and (b) c-axes, respectively. The solid (dashed) line 
corresponds to the low-temperature specific heat calculated in 
the framework of the sine-Gordon (boson) model and given 
in Eq. QJj] (Eq. O . Dotted line plotted in (a) corresponds to 
aT 3 lattice contribution to the specific heat. 



Neel temperature increases with increasing field strength, 
as shown in Fig. Hfb). 

We analyze the specific heat for H > H c to deter- 
mine the field dependence of the spin gap. The low- 
temperature specific heat is first compared to a simple 
model (which we refer to as the boson model) of n species 
of one-dimensional non-interacting gapped bosons with a 
dispersion relation 



hu>(q) = y/A 2 + [v{q 



9o J 



(5) 



Here A is the spin gap, v is the spin-wave velocity, which 
depends on both the strength and orientation of the mag- 
netic field, q is the chain wave- vector q = irh and qo = tt. 
The specific heat of a one-dimensional boson gas in this 
boson model is given by 15 



Cmag (T) 



hR 



A 



V2^ \k B T 



exp(-A/£; B T) . (6) 



Terms proportional to -yj- and to T 3 were included in the 
fit to take into account the small nuclear spin and lattice 
contributions, respectively. The total function fitted to 
the specific heat was 



C{T, H) = C mag + aT 3 + b{H/Tf 



All data for T < A and H > H c were fit simulta- 
neously to Eq. yielding a = 0.027(1) J/molK 4 and 
b = 8(2) • 10~ 6 JmolK/T 2 . The average spin-wave veloc- 
ity is found to be v/h = 0.65(1) and 0.58(2) meV/mode 



4 



for fields along the a- and c-axis, respectively. From 
our previous measurements of the field-induced incom- 
mensurate mode^ we found that the spin-wave velocity 
is v = 1.84(4) meV. This yields n = 2.8 for the number 
of gapped low-energy modes which contribute to the spe- 
cific heat, in excellent agreement with the expectation of 
one longitudinal and two transverse modes. The result- 
ing field and orientation dependent spin gaps are shown 
in Fig. [4] and the calculated C(T, H) for the boson model 
are plotted as dashed lines in Figs. [5] and [31 
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FIG. 4: Field dependence of the energy gaps for the boson and 
sine-Gordon models, derived from fits to specific heat data as 
shown in Fig. [3] The error bars are smaller than the symbol 
size. 



where Mo = M. We use this expression for an additional 
comparison to the experimentally-observed temperature 
and field dependence of the specific heat C(T,H) for T < 
Mi and H > H c , in order to determine the soliton mass 
as a function of field strength and orientation. These 
results are shown as solid lines in Figs. and O The 
adjusted soliton and breathers masses are shown in Fig.[2J 
in comparison with the results from the boson model. 

The adjusted boson mass is larger than the breather 
mass, but smaller than the soliton mass, as expected from 
a model which is only sensitive to an average gap energy. 
The soliton energy is different for fields applied along 
the a- and c-axis as a result of different staggered fields 
arising from the staggered g- factor and DM interactions. 

The field dependence of the gap for two different field 
directions allows an estimate of the strength of g s and 
D/J. The best estimate is obtained for high fields where 
the effect of the low-field ordered phase is small. Assum- 
ing H st = cH and fitting Eq.[5]to the field dependence of 
the gap, we obtain c = 0.045(2) for fields along the a-axis 
and c = 0.090(4) for fields along the c-axis. Comparing 
these results with Eq. Q] we obtain 

g s - g aa D/2J = 0.045 

g s + g cc D/2J = 0.09. (11) 

With g aa = 2.28 and g cc = 2.12 (Ref. this leads 
to g s = 0.068(3) and D/J = 0.0102(5), establishing the 
presence of both sizeable DM interactions and a stag- 
gered g factor in CDC. 



We also analyze the specific heat in the framework of 
the sine-Gordon model, to which the Hamiltonian can 
be mapped in the long wave-length limit^. The sine- 
Gordon model is exactly solvable, with solutions that 
include massive soliton and breather excitations. The 
soliton mass, M, is given b y 16 i 17 



J 



M 

\ (2^-/3^)74^/2 



£i)l/4T-(l+£)/2_ 



(8) 



Here £ = f3 2 /(8n - f3 2 ) -> 1/3 for /j, H 0, B = 
0.422169, and /J,oH st is the effective staggered field. The 
mass of the breathers, M n , is given by 



M n = 2M sin(n7r£/2) 



(9) 



where 1 < n < [l/£] with £ depending entirely on the ap- 
plied field^. The temperature dependence of the specific 
heat due to breathers and solitons is given by 
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FIG. 5: (h,0,0) scans through Q = (1,0,0) magnetic Bragg 
peak in CDC at T = 0.25K, n H = 7.25 T||b, T = 0.25K 
and T = 2K at zero field. The bar in the figure represents the 
calculated peak width of Q = (1,0,0) based on the known 
instrumental resolution. The peak at h = 1.09 appears to be 
field and temperature independent indicating that it does not 
arise from magnetic scattering. 
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V. NEUTRON DIFFRACTION 

Neutron diffraction measurements were performed in 
the (h,0,l) and (h,k, 0) reciprocal lattice planes. Fig. [5] 
shows /i-scans through the (1, 0, 0) reflection for different 
temperatures and fields. At zero field and T = 0.3 K a 
resolution- limited peak is observed at h = 1. The peak 
is absent for temperatures higher than T)v = 0.93 K. 
No incommensurate order was observed in any of our 
experiments. The reflections we found can be indexed 
by (2 • n,m, 0) in the (h,k,0) plane, whereas magnetic 
diffraction was found only at the (1, 0, 0) and the (3, 0, 0) 
reciprocal lattice points in the (h,0,l) plane for |Q| < 
2.6A- 1 . 
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FIG. 6: Measured field dependence of the Q = (1,0,0) mag- 
netic peak intensity at T = 0.25 K with the field applied along 
the c-axis and at T = 0.25K and at T = 1 K with the field 
applied along the b-axis. The solid lines represent the fits 
explained in Sect. IV CI 



A. HT-Phase Diagram 

The phase diagram of CDC was further explored with 
neutron diffraction for fields up to 7.25 T, applied along 
the b and c-axis. Figure [S] summarizes the magnetic 
field dependence of elastic scattering measurements at 
the (1,0,0) wave-vector. Applying a field along the b- 
axis increases the intensity of the (1,0,0) reflection and 
we associate this with an increase in the staggered mag- 
netization. In contrast, a field along the c-axis decreases 
the intensity of the (1,0,0) Bragg peak, first in a sharp 
drop at 0.3 T suggestive of a spin-flop transition, and 
then in a continuous decrease towards a second phase 
transition at 3.9 T. The transition at 3.9 T is evidence 
that exchange-induced long-range order is suppressed by 
the application of fields along the c-axis. The value of 
the critical field, 3.9 T, is in agreement with specific heat 
measurements. 



FIG. 7: Measured temperature dependence of the Q = 
(1,0,0) magnetic peak intensity for zero field, fioH = 7.25 T 
applied parallel to the b-axis and /j,qH — 3 T applied parallel 
to the c-axis. The solid lines represent the fits explained in 
Sect. |VCl 



Fig. [7] shows the temperature dependence of the 
(1, 0, 0) magnetic Bragg peak for fields applied along the 
b- and c-axis. For fields along the b-axis, the Bragg peak 
has higher intensity at low T and the enhanced intensity 
survives to higher temperatures than in the zero field 
measurement. For fields along the c-axis, the intensity 
is reduced and the transition temperature is lower than 
at zero field. The phase transition remains continuous, 
regardless of the direction of the field. 
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FIG. 8: Field-temperature phase diagram of CDC determined 
from heat capacity measurements (filled symbols) and from 
neutron scattering (open symbols). The squares indicate the 
spin-flop transition observed for a field along the c-axis. The 
solid and dashed lines are a guide to the eye. 



The phase diagram in Fig. [8] shows a synopsis of the 
results obtained from specific heat and neutron diffrac- 
tion measurements. CDC adopts long-range magnetic 
order below Tjv. This order can be strengthened by the 
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application of a magnetic field along the b-axis, which en- 
hances the ordered moment and increases the transition 
temperature to long-range order. In contrast, applying 
fields along the a- and c-axes competes with the zero-field 
long-range order, decreases the critical temperature lead- 
ing eventually to the collapse of the low-field long-range 
order at 5.9 and 3.9 T respectively, as extrapolated to 
zero temperature. 



B. Magnetic Structures 



ble U and Fig. [TOb-b) with the experiment is obtained 
for a magnetic structure belonging to representation T 2 
as defined in Appendix IB1 with \ 2 = 3.1 and R — 0.43. 
The magnetic structure is collinear with a moment of 
0.44(5)/is along the c-axis as shown in Fig. [9ja). The 
alignment of the moment along the c-axis suggests an 
easy-axis anisotropy along this direction, which is not in- 
cluded in the zero-field Hamiltonian defined by Eq. [T] 
The ordered magnetic moment is significantly less than 
the free-ion moment, as expected for a system of weakly- 
coupled S = A chains close to a quantum critical point 19 . 



To understand the field-induced phase transitions, it 
is important to determine the symmetry of the mag- 
netic structures below and above the transitions. The 
intensities of ~ 40 magnetic Bragg peaks were measured 
at 0.1 K at zero field, at fi H = 0.5 T and at 10 T 
applied along the c-axis. The ordering wave-vector for 
both structures is k = (0,0,0). We use group theory 
to identify the spin configurations consistent with the 
wave-vector k for the given crystal symmetry, and com- 
pare thestructure factor of possible spin configurations 
with the observed magnetic Bragg peak intensities. The 
group theoretical analysis which leads to the determina- 
tion of the symmetry-allowed basis vectors is presented 
in Appendix [Bl 






FIG. 9: Magnetic structures at zero field (a), 0.75 T (b) and 
10 T (c), with the field H applied along the c-axis. The solid 
arrows indicate the spin directions, + and - indicate spins 
pointing into opposite directions along the b-axis. The Cu 2+ 
positions are defined in Appendix [A] Because the ordering 
wave- vector is k = (0,0,0), the spin arrangement in other 
unit cells is identical, (d) Directions of the staggered fields 
H st upon application of an external field H along the c-axis. 
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TABLE I: Measured magnetic Bragg peak intensities at 
HoH = T, compared with the best fit with % 2 = 3.1 and 
R — 0.43. The found magnetic structure belongs to T 2 and 
the moment is aligned along the z-axis. The next best fit 
yields a magnetic structure that belongs to T 5 with the mag- 
netic moments aligned along the y-axis, with \ 2 = 5.2 and 
R = 0.67. 

At ixqH = 0.75 T after the system undergoes a spin- 
flop transition, the magnetic order can be described by 
the representation T 5 with x 2 — 1-75 and R = 0.49 ( 
Table |TT] and Fig. [TUb ) . Slightly lower values for \ 2 are 
obtained if a small antiferromagnetic moment belonging 
to r 8 is also allowed for. As shown in Fig. [DJb), the 
spin structure is collinear with a magnetic moment of 
0.3(1)/ib per Cu 2+ along the b-axis (and perpendicular 
to the field as expected for a spin- flop transition). 

At fioH = 10 T, the magnetic order can be described 
by representation T 8 with \ 2 = 2.41 and R = 0.21 ( Ta- 
ble IIIII and Fig. flU)) . The spin structure is a collinear 
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TABLE II: Measured magnetic Bragg peak intensities at 
fioH = 0.75 T, compared with the best fit with x 2 = 1-75 
and R — 0.49. The found magnetic structure belongs to T 5 
and the moment is aligned along the y-axis. The next best 
fit yields a magnetic structure that belongs to T 4 with the 
magnetic moments aligned along the x-axis, with \ 2 = 5.3 
and R = 0.96. 



magnetic structure where the magnetic moments point 
along the a-axis, perpendicular to the external field and 
along the staggered fields in the material, as shown in 
Fig. [He). The antifcrromagnctically ordered moment 
along the a-axis is 0.49(5)/zb per Cu 2+ . Our numerical 
calculations- using a the Lanczos method for finite chains 
of 24 sites yield a ordered staggered moment 0.55/is per 
site, in good agreement with the experiment. 

The magnetic structures reveal why fields along the a- 
and c-axis lead to the collapse of long-range order, while 
a field along the b-axis strengthens magnetic order. A 
field along the c-axis leads to staggered fields due to the 
alternating g tensor and DM interactions, as illustrated 
in Fig. Hfd). These staggered fields compete with inter- 
chain interactions which favor a different antiferromag- 
netic arrangement of spins. For example, inter-chain ex- 
change favors AF alignment of the neighboring magnetic 
moments on sites 1 and 4, but the staggered fields favor 
ferromagnetic alignment along the a-axis. This must lead 
to an increase in quantum fluctuations which leads to a 
collapse of magnetic order, identifying the transition as 
a quantum phase transition. This situation was recently 
analyzed in great detail by Sato and Oshikawa^. 

Application of a field along the b-axis does not induce 
staggered fields, and so merely quenches quantum fluctu- 
ations for small field strengths, thereby leading to an en- 
hancement of the antiferromagnetically ordered moment 
as observed in the experiment. Application of a field 
along the a-axis should not lead to a spin-flop transition, 
because the moments in the zero-field structure are al- 
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TABLE III: Measured magnetic Bragg peak intensities at 
HoH = 10 T, compared with the best fit with \ 2 = 2.4 and 
R = 0.21. The found magnetic structure belongs to T 8 and 
the moment is aligned along the x-axis. The next best fit 
yields a magnetic structure that belongs to T 6 with the mag- 
netic moments aligned along the z-axis, with \ 2 = 8.2 and 
R = 0.34. 

ready perpendicular to that field direction. At higher 
fields, however, the competition between magnetic order 
arising from either spin exchange or staggered fields be- 
comes important, leading eventually to the phase transi- 
tion at fioH = 5.7 T into a gapped phase that is described 
by the sine-Gordon model. 

C. Critical Exponents 

The neutron diffraction measurements allow determi- 
nation of the order parameter critical exponent (3 for 
the phase transition to long-range order. Fig. [TT1 shows 
the magnetic intensity of the (1,0,0) reflection close to 
the critical region. The data were corrected for a non- 
magnetic, constant background. 

The temperature scans are shown as a function of re- 
duced temperature (1 — T/T c ) for zero field, (j,qH = 
7.25 T applied along the b-axis and n^H = 3 T applied 
along the c-axis. Power law fits were performed using 

1(H) oc \T - T N \ 20 ■ (12) 
It was found that (3 = 0.30(1) for zero field, (3 = 0.35(1) 
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FIG. 10: (a-d) Observed magnetic intensities plotted as a 
lunction of calculated intensities, showing the quality of our 
fits, (a-b) compares the best fit for ^qH = T with the 
second best fit, while (c-d) show the best fits for fj,o H — 0.75 T 
and fj,oH = 10 T, respectively. 



for fi H = 7.25 T and (3 — 0.36(1) for ^ H = 3 T. 
were performed over a reduced temperature of 0.2. 
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FIG. 11: Logarithmic plot of the background corrected Q = 
(1,0,0) magnetic intensity plotted as a function of reduced 
temperature (open symbols) or versus reduced field (filled 
symbols). Dashed lines corresponds to the power law fits de- 
scribed in the text. For the temperature scans, fioH = 7.25 T 
was applied along the b-axis, and fioH = 3 T was applied 
along the c-axis. For the field scans at constant temperature, 
the field was applied along the b-axis for T = 1 K and along 
the c-axis for T = 0.25 K. 



Figure QT] also shows two field scans: one scan with the 
field applied along the c-axis and performed at 0.25 K, 
and a scan with the field applied along the b-axis, per- 
formed at 1 K. The latter transition is unlike all other 



transitions discussed here, because here the field induces 
an ordered state starting from a paramagnetic state. 
The field scans are shown as a function of reduced field 
|1 — H/H c \. Power law fits were made using 



J(T)« \H-H C \ 2 ? 



(13) 



in order to determine the critical exponents and it was 
found that j3 = 0.400(1) for T = 0.25 K and (3 = 0.481(3) 
for T = 1 K. Fits were performed over the range 0.2 in 
reduced field. 

For the temperature driven transition, the critical ex- 
ponents (3 fall in the range 0.30 ~ 0.36 depending on the 
direction of the applied field, consistent with the critical 
exponents of three-dimensional Ising universality class^i 
to which CDC belongs at zero field. For the field-induced 
transitions, the Ising character of the spin symmetry is 
enhanced due to the presence of a uniform and perpendic- 
ular staggered field. The critical exponent for the field- 
induced transitions, (3 = 0.400(1) and (3 = 0.481(3), are 
however higher than the critical exponent (3 = 0.325 of 
the three-dimensional Ising universality class. 

For an Ising quantum phase transition (T = 0), the 
effective dimension of the quantum critical point is D = 
d + z where d — 3 is the spatial dimension and z = 1 is 
the dynamical exponent. Since D — 4 is the upper crit- 
ical dimension, the quantum critical point is Gaussian 
(with logarithmic corrections) and (3 = 0.5. Therefore, 
by measuring the field dependence of the order param- 
eter at different temperatures, we expect to observe a 
crossover between the classical (finite T and (3 — 0.325) 
and the quantum phase transitions (T = and (3 = 0.5). 
This seems to be the natural explanation for the observed 
values of f3 close to 0.5. A similar increase in the appar- 
ent f3 critical exponent was observed in the singlet ground 
state system PHCC^ 



VI. INELASTIC NEUTRON SCATTERING 

A. Exchange Interactions 

To characterize the Hamiltonian and accurately deter- 
mine the strength of the spin interaction along the chain, 
we performed an energy scan for wave-vector transfer 
Q • a = 7r along the chain. At this wave-vector the two- 
spinon spectrum characteristic of the AF S — \ chain 
should be concentrated over a small energy window and 
form a well-defined peak located at hco = ttJ/2. Fig. [T2l 
shows that this is indeed the case, and that for CDC this 
well-defined excitation is located at fiuj — 2.43(2) meV. 
A fit to the exact two-spinon cross-section^ convolved 
with the experimental resolution yields an excellent de- 
scription of the observed spectrum, as shown in Fig. [T^] 
for a value of intra-chain exchange J = 1.46(1) meV. 
This value is very similar to that obtained from the 
temperature dependence of the magnetic susceptibility, 
which was 1.43 meVJ^ 
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FIG. 12: Neutron scattering intensity as a function of en- 
ergy transfer for zero-field for h = 0.5 measured using SPINS. 
The measurement was performed with 9 analyzer blade chan- 
nels set to reflect Ef = 5 meV to the center of the detec- 
tor. A Be filter rejected neutrons with energies higher than 
5 meV from the detection system. The measurements were 
performed with the strongly dispersive direction along the 
scattered neutron direction, so that the measurements inte- 
grate neutron scattering over wave-vectors along a weakly- 
dispersive direction. The solid line represents the exact two- 
spinon cross-section^ for J = 1.46(1) meV, convolved with 
the experimental resolution. 



B. Ground state energy at zero field and 

HoH = 11 T 

Inelastic neutron scattering not only allows determi- 
nation of exchange constants, but also a measurement 
of the ground state energy. The expectation value of 
the Hamiltonian, < TL >, is proportional to the first 
moment of the dynamic structure factor, defined as 
(lu)q = H 2 J du>u)S(Q,u}). Here the normalization of 
S(Q, uS) is chosen such that the total moment sum rule 
reads: h J 5(Q, uj)dui = S(S + l)/3. So if the dynamic 
structure factor, S(Q,u>), is measured in absolute units, 
neutron scattering yields a measurement of < TL > as a 
function of field or temperature. 

For a Heisenberg AF chain, the wave-vector depen- 
dence of (uj) h is given by 

{Cj)h = -\ < H> (1 - cos(Trft)) . (14) 

We use this relation to describe the measured first mo- 
ment, (uj)h- Figure [l~3l shows that at zero field the ex- 
perimental (uj)h follows this expression with < TL >~ 
—0.4(1) J - in excellent agreement with Bethe's result of 
< TL >~ -0.44J (Ref. |H). {Cb) h was obtained by fitting 
the exact two-spinon cross-section 5 , 2 Sp (Q,^) as calcu- 
lated by Bougourzi et al. 23 to the data and computing 
(ui)h from the analytical expression for 5*2 S p(Q,w). 

Fig. [1~3"1 shows the first moment (u>)h at /i H = 11 T in- 
cluding contributions only from sharp modes in the spec- 



1.00 




h (r.l.u.) 

FIG. 13: First moment {u))h ~ h 2 J dujujl(h,ui) of the well- 
defined excitations only as a function of wave-vector h along 
the chain, for zero field and hqH=11 T. The raw scattering 
data were corrected for a time-independent background mea- 
sured at negative energy transfer, for monitor efficiency, and 
for the Cu 2+ magnetic form factor, folded into the first Bril- 
louin zone, and put onto an absolute scale using the elastic 
incoherent scattering from CDC. 



trum, thus excluding any contributions from continuum 
scattering. It was calculated by fitting a Gaussian to the 
experimentally observed peak in the energy spectra and 
by integrating the area under the Gaussian curve. (u>)h 
at HqH = 11 T is generally about 3 times smaller than 
(u))h at zero field. This is evidence for the presence of 
substantial magnetic continuum scattering because both 
numerical and analytical calculations show that the first 
moments in an applied magnetic field should be larger 
than this measured value. This is based upon (A) nu- 
merical calculations^ using the Lanczos method for finite 
chains of length between 12 and 24 spins which show that 
the ground state energy at fioH = 11 T for H st = 0.075-ff 
is < TL > = —0.34 J, and (B) analytical calculations of the 
field dependence of the first moment which show that 
Eq. [H] acquires merely a small constant value with the 
application of uniform and staggered fields. Points (A) 
and (B) thus show that theoretically, only a slight depres- 
sion of the first moment is expected upon application of 
magnetic fields. 

Because our data presented in Fig. ll3l include only res- 
onant modes this result indicates that there are substan- 
tial contributions to the first moment from continuum 
scattering - perhaps not a surprise given that continuum 
scattering is the sole source of the first moment at zero 
field. 



C. Comparison with the chains-MF model 

The magnetic order of weakly-coupled chains can be 
described by a chain mean-field theory developed by 
Schulz —. It gives the following relations between ordered 
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magnetic moment Too , Tjv and the average of inter-chain 
interactions \J'\: 



1.28y/]n{5.8J/T N ) 

m = l.017y/\J'\/J. (16) 

The first relation allows an estimate for |J'|. In our 
case, where Tn = 0.93 K and J = 1.5 meV, this gives 
| J'| = 0.03 meV. The real value for |J'| may be some- 
what higher because quantum fluctuations are not fully 
accounted for in the chain mean-field theory. 

The second relation predicts the ordered moment for 
a Heisenberg chain magnet too = 0.14 [1b for \J'\ = 
0.03 meV, which is considerably lower than the observed 
ordered moment even if one considers a higher value 
for | J' | due to quantum fluctuations. The large dis- 
crepancy may be due to the Ising-type spin anisotropy 
which fixes the zero-field ordered moment along the c- 
axis, thereby quenching quantum fluctuations and en- 
hancing long-range order at low temperatures. Indeed, 
this term, no matter how small, induces long-range order 
even in the purely ID model at T = K. 



was characterized by an apparent critical exponent /?, 
which is higher than expected for a classical phase tran- 
sition in a 3D antiferromagnet. The interpretation of this 
result is that the experiment probes a cross over regime 
affected by the T = quantum critical transition, which 
is above its upper critical dimension and therefore mean 
field like. 
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APPENDIX A: DEFINITION OF THE 
MAGNETIC SUB-LATTICE IN CDC 



VII. CONCLUSIONS 

CDC is a rich model system in which to study quan- 
tum magnetism because an applied magnetic field leads 
first to a first-order phase transition and then to a novel 
high-field phase with soliton excitations via a quantum 
critical phase transition. In this paper, we character- 
ized the material and its proximity to quantum critical- 
ity using specific heat and neutron scattering measure- 
ments. We determined the phase diagram as a function 
of temperature and field applied along all three crystal- 
lographic directions. The Bragg peak intensities and the 
temperature dependence of the specific heat demonstrate 
that when applying a field along the a- and c-axes, the 
spin system undergoes a quantum phase transition due 
to competing magnetic order parameters. Beyond the 
critical field, specific heat measurements show activated 
behavior typical for a gapped excitation spectrum. The- 
oretical analysis indicates that the spin gap results from 
spinon binding due to a staggered gyromagnetic factor 
and DM interactions. The gap energy was extracted 
from specific heat data using both a non-interacting bo- 
son model and an expression derived from a mapping of 
the Hamiltonian to the sine-Gordon model. We deter- 
mined the strength of the staggered gyromagnetic fac- 
tor and the DM interactions by analyzing the different 
field dependence of the gap for fields along the a- and 
c-axes. A magnetic field along the c-axis leads to two 
phase transitions: a first-order spin-flop transition and a 
continuous phase transition at higher field to magnetic 
order described by only one irreducible representation. 
The continuous field driven quantum phase transition 



The unit cell contains four Cu 2+ sites, which occupy 
the 4c position in the Pnma space group and are given 
by: 

T\ = (0.3209,0.25,0.3614) 
r 2 = (0.8209,0.25,0.1386) 
r 3 = (0.1791,0.75,0.8614) 
r 4 = (0.6791, 0.75, 0.6386) (Al) 

This numbering of the Cu 2+ positions is used throughout 
the paper. 

APPENDIX B: MAGNETIC GROUP THEORY 
ANALYSIS 

Magnetic structures were inferred from diffraction data 
by considering only spin structures "allowed" by the 
space group and the given ordering wave-vector k — 
(0,0,0). Landau theory of continuous phase transitions 
implies that a spin structure transforms according to an 
irreducible representation of the little group Gk of sym- 
metry operations that leave the wave- vectors k invariant. 
The eigenvectors <fi x of the A-th irreducible representa- 
tions r A were determined using the projector method. 25 
They are given by 

^ A = E^(.9)9W' ( B1 ) 

9 

where g is an element of the little group and <p is any vec- 
tor of the order parameter space. x A (<?) is the character 
of symmetry element g in representation T A . 
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TABLE IV: Basis vectors of irreducible representations for the 
commensurate phase described with k = (0, 0, 0) for the for 
Cu +2 sites. The components of the vector correspond to the 
spin component on the Cu sites in order given in Table |A"T1 

We start by considering the symmetry elements of the 
Pmna space group of CDC: 

{1, T, 

2 £ ) 2y j 2 2 , Tfl X y , VTl xz , Tflyz 

}■ (B2) 

Here 1 is the identity operator, 1 is inversion at the origin, 
2 Q denotes a 180° screw axis along the crystallographic 



direction a — a, b or c (180° rotation operation followed 
by a translation) . m, a p is a glide plane containing axes a 
and (3. The group is nonsymmorphic because the group 
elements {-R|a} consist of an operation R followed by a 
translation a equal to half a direct lattice vector. 

The ordering wave- vector k = (0, 0, 0) is invariant un- 
der all operations of the space group, so the little group 
Gk equals the space group. The group consists of 8 dif- 
ferent classes and therefore has 8 irreducible represen- 
tations. The irreducible representations and their basis 
vectors are given in Table [TVl The magnetic representa- 
tion can be written as 

r k = r 1 e 2 r 2 ® r 3 ® 2 r 4 
®r 5 ®2r 6 ®r 7 ®2r 8 . (B3) 



The low-temperature magnetic structure was inferred 
from the integrated intensities of rocking scans through 
magnetic Bragg peaks in the (h, k,0) plane. An abso- 
lute scale for the integrated intensities of these reflec- 
tions was obtained from measurements of nuclear Bragg 
peaks. The magnetic structure factors squared were ob- 
tained after correcting for resolution effects, the sample 
mosaic and the magnetic form factor of the Cu 2+ ions. 
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